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I. INTRODUCTION 



For a closed system of two non-relativistic point masses interacting via a central force the 
angular momentum L of the relative motion is conserved and the motion is always in the 
plane perpendicular to L. If the force is derived from the Newton or Coulomb potential, 
there is an additional conserved quantity: the Laplace- Runge-Lenz 1 vector ffl. For the 
hydrogen atom this vector has the form 

1 e 2 
C = — p x L r , L = r x p , 

m r 

where m denotes the reduced mass of the proton-electron system. The Laplace-Runge-Lenz 
vector is perpendicular to L and hence is a vector in the plane of the orbit. It points in the 
direction of the semi-major axis. 

Quantum mechanically, one defines the hermitian Laplace-Runge-Lenz vector 

1 e 2 
C = — (pxL-Lxp) r. (1) 

2m r v 1 

By exploiting the existence of this conserved vector operator, Pauli calculated the spec- 
trum of the hydrogen atom by purely algebraic means [|3], |4[]. He noticed that the angular 
momentum L together with the vector 



which is well-defined and hermitian on bound states with negative energies, generate a 
hidden 5*0(4) symmetry algebra, 

[L a , Lb] = ihe abc L c , [L a , K b \ = ihe abc K c , [K a , K b \ = ihe abc L c , 

and that the Hamiltonian can be expressed in terms of K 2 + L 2 , one of the two second-order 
Casimir operators of this algebra, as follows 

me 4 1 

~~K 2 + L 2 + fr 2 ' ( 2 ) 
One further observes that the other Casimir operator K • L vanishes and arrives at the 
bound state energies by purely group theoretical methods. The existence of the conserved 



1 A more suitable name for this constant of motion would be Hermann-Bcrnoulli-Laplace vector, see 
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vector K also explains the accidental degeneracy of the hydrogen spectrum. Only much 
later the scattering amplitude of the hydrogen atom has been calculated with the help of 
the Laplace- Runge-Lenz vector ||. 

In this paper we shall generalize these results in two directions: first to the hydrogen atom 
in arbitrary dimensions 2 || and second to the corresponding supersymmetric extensions. 
In the following section we prove that the Schrodinger Hamiltonian in d dimensions with 
1/r potential admits a generalization of the Laplace- Runge-Lenz vector. Together with the 
generators of the rotation group SO(d) this vector generates the dynamical symmetry group 
SO(d+ 1). This hidden symmetry allows then for a purely algebraic solution of the hydrogen 
atom in arbitrary dimensions, very much as in three dimensions. 

In section [HI] we summarize the extensions of <i-dimensional Schrodinger Hamiltonians to 
models with M = 2 supersymmetry. The corresponding Hamiltonians may be written 
2 d -dimensional matrix Schrodinger operators. For a particular choice of the superpotential 
we obtain the supersymmetric extension of the hydrogen atom. For instance in d = 3 (and 
in a suitable basis) we find the Hamiltonian 

/l \ 

Mi 

M 2 



with 3x3 matrices 



Mi ab = X a X b , M 2a b = Sab ~ {-) a+b X a X b 



In section [TV] we construct the supersymmetric extensions of the angular momentum and the 
Laplace-Runge-Lenz vector. Similarly as for the purely bosonic system, together they form a 
dynamical SO(d+l) symmetry algebra. This symmetry is exploited in the following section 
to obtain the discrete eigenvalues and their degeneracies. In section [VI| we characterize all 
bound state wave functions. In the last section we illustrate our general results by analyzing 
in detail the supersymmetric hydrogen atoms in two, three and four dimensions. 
The appendices contain the necessary group theoretical tools needed in the main body of 
the paper. 



2 When speaking of the d— dimensional hydrogen atom, we always mean the 1/r-potential, although this 
potential permits the application of Gauss' law in three dimensions only. 
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II. THE COULOMB PROBLEM AND ITS SYMMETRIES IN d DIMENSIONS 



We generalize the Coulomb problem to arbitrary dimensions by keeping the 1/r-potential, 
although this potential solves the Poisson equation in three dimensions only With this 
assumption the hydrogen atom in d dimensions is governed by the Schrodinger equation 

m ' = {-&n A ~t)* = ' GO 

It is convenient to measure distances in units of the Compton wavelength X c = h/mc. With 
respect to these dimensionless coordinates equation @ takes the simpler form 

Ht/j = E^, H = p 2 -^, p a = -d a , a = l,...,d, (4) 

r 1 v ' 

where rj is twice the fine structure constant a and the dimensionless energy E is measured 

in units of mc 2 /2. The central force is attractive for positive rj. 

The hermitian generators L a b = x a p b — x b p a of the rotation group in d dimensions satisfy 
the familiar so(d) commutation relations 

[L a b, L cd ] = i(S ac L bd + SbdL ac — S a dL bc — 5b c L ad ) , (5) 
where indices run from 1 to d. It is not very difficult to guess the generalization of the 
Laplace- Runge-Lenz vector (|l|) in d dimensions ||, 

TJX a 

C a = L ab p h + p b L ab — . (e) 

Indeed, these operators commute with the Hamiltonian (Q). They form a SO (rf)-vector, 

[L a b, C c ] = i(5 ac C b — 5 bc C a ) , 
and the commutator of C a and C b is proportional to the angular momentum: 

[C a ,C b } = -A\L ab H . (7) 

Now we may proceed as we did in three dimensions and define on the negative energy (E < 0) 
subspace of the Hilbert space H = L 2 (]R d ) the hermitian operators 

K a = ^-^jj c a , with [K a , K b ] = iLab ■ (8) 
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The operators {L ab , K a } form a closed symmetry algebra of dimension (d + l)d/2. More 
explicitly, they can be combined to form generators Lab of the orthogonal group 3 SO(d + 1) 



J AB 



L a b 


K a 


-K b 






(9) 



which implies that the Lab obey the commutation relations (|5|) with indices running from 
1 to d + 1. We can verify a relation similar to (0) by solving 

C a C a = —4:K a K a H = V 2 + (2L ab L ab + (d - lf)H (10) 
for H . We obtain the Hamiltonian 

it - 2 *7 _ V 2 

H ~ P "-"(d-l)* + 4C (2) (") 
in terms of the second-order Casimir operator of the dynamical symmetry group SO(d+ 1), 



C(2) — -LabLab — -L a bL a b + K a K a . 

This is the generalization of @ we have been looking for. It remains to characterize those 
irreducible representations of SO(d+ 1) which are realized in the Hilbert space L 2 (M <i ). In 
three dimensions the allowed representations are fixed by the condition K • L = on the 
second Casimir operator of SO (A). We expect n — 1 conditions on the n Casimir operators 
of the dynamical symmetry group SO(2n+l) in d = 2n dimensions and n conditions on the 
n+1 Casimir operators of SO (2n+2) in d = 2n+l dimensions. In the following we treat the 
even- and odd-dimensional cases separately. 

Even-dimensional spaces: An irreducible representation is uniquely characterized by 
its highest weight state. By definition, this state is annihilated by all raising operators 
belonging to the simple roots. To characterize these states one conveniently chooses complex 
coordinates zi,...,z n in R 2n such that the raising operators of the dynamical symmetry 
group with generators Jab i n © have the simple form (cf. appendix 0) 

Ei = -(zid i+1 - z i+1 di) , i = 1, . . . , n - 1 , (12) 

E n = -^(Cd-i + iC d ) = -2z n A + (2rd r + d - l)d n - ^ . (13) 



3 For scattering states (E > 0) a similar redefinition leads to generators of the Lorentz group SO(d, 1). 
Here we are interested in bound states and will not further discuss this possibility. 
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In the formula for E n we actually should have used the operators Kd-i and Kd in @, since 
they appear as components of Jab- But since we are only interested in highest weight states 
which are annihilated by E n we may take the operators C^-i and instead. Also note 
that all simple roots of SO(d) are positive roots of SO(d + l) so that all highest weight 
states of SO(d + l) are automatically highest weight states of SO(d) C SO(d+l). Now it 
is not difficult to see that a regular wave function which is annihilated by all simple raising 
operators of SO(d), that is by the n—1 operators in (|l^) and by the operator z n _id n — z n d n -i 
(see appendix |A 1|) , must have the form (cf. ( |A10| )) 

* = f(r)z{ . 

It is a highest weight state of SO(d + 1), if in addition it is annihilated by E n in (|T3|): 

E n V = - |(d - 1 + 2£)^log/ + ^|^ = 0. (14) 
Hence the highest weight state reads 

The constant r\ must be positive for bound states to exist. ^ is of course an eigenfunction 
of all n Cartan generators Hi = Zidi — Zidi of SO(2n + 1) with eigenvalues (£, 0, . . . , 0). That 
shows that only the symmetric multiplets of the dynamical symmetry group appear 4 . From 
we take the values of the second-order Casimir operator for symmetric multiplets 



C (2) =£(£ + d-l) , £ = 0,1,2,... (16) 
and their dimensions 

+ d\ f£ + d-2' s 



dim Vf 



- 2 



(17) 



The dimV^ states of the symmetric representation are obtained by acting repeatedly with 
the lowering operators 

El <n = l{z i+1 d i -z i d i+1 ) and E\ = -2z n A + d n (2rd r - 1) - ^ 
l r 

on the state (O). This way one obtains all bound if-eigenstates with the same energy. 



This corresponds to the extension of Fock's method to d dimensions, cf. j|, 0, 



Odd-dimensional spaces: For d = 2n+l the rank of the dynamical symmetry group 
SO(d + 1) exceeds the rank of the rotation group SO(d) by one. We combine the first 2n 
coordinates to n complex ones and take coordinates Z\, . . . , z n , x d in R d , see appendix |A2 . 
As Cartan generators we choose 

Hi = Zidi - zA , i = 1, . . . , n and H„ +1 = C d . (18) 
The raising operators are the n — 1 operators E{ in (|1^) supplemented by 

E n = \{z n d Xd - x d d Zn ) and E n+1 = ^=(C d _ 2 + iQ-i) • 
i v 2 

The last raising operator E n+ i coincides with E n in (|T3|). A regular wave function is anni- 
hilated by the first n raising operators only if it has the form ^ = f(r)z{. The requirement 
that it is annihilated by the last raising operator E n+ i again leads to equation fll4j) and 
hence to the solution \1/ in (|15|). To determine the multiplets with this highest weight state 
we need to calculate the highest weight vector, that is the value of the Cartan generators 
on ty. Clearly, 

H x m = and = for i = 2, . . . , n . 

The last Cartan generator in (|18D has the explicit form 

H n+1 = ~2x d A + (2rd r + d- l)d 3 

i 

and we find 

#n+i* = - { (d - 1 + 2£) A log / + V j ^ = , 

after using / = e~ ler . Hence, on any highest weight state the operators H 2 , . . . , H n+1 vanish 
and again we find the completely symmetric representations of the dynamical symmetry 
group SO(d+ 1). The eigenvalues of the second-order Casimir operator and the dimension 
of the representations are given by the same formulae ( |16| , [iT| ) as for the hydrogen atom in 
even dimensions. 

Since rj is twice the fine structure constant a = e 2 /he and E is measured in units of mc 2 /2 
the formula (|ll|) yields the well-known |§ bound state energies in three dimensions 



rjx d 



The degeneracy of E n is the dimension n 2 of the corresponding symmetric representation of 
SO (A). All n 2 states with the same energy E n are gotten by acting with the two lowering 
operators on the highest weight state 



Xi + ix 2 ) 



a 

In = ~ 

n 



In d dimensions the corresponding formulae read 



Ep 



me 



+ 



SSf(x) 



Xi + IX2) 



a 



£+(d-l)/2 



2K 2 V ' 2 
and Eg has degeneracy dim Vg in (|I7|). 

The appearance of the accidental degeneracy - phrased in the language of representation 
theory - corresponds to the following branching rule: the completely symmetric representa- 
tions of the dynamical symmetry group SO(d+l) branches into those completely symmetric 
representations of the rotation group SO(d) with equal or shorter Young diagrams, 



1 


2 






t 



1 


2 






I 



SO(d) 



(19) 



SO(d+l) 

all of them possessing the same energy. The energy, its degeneracy and the bound state 
wave functions are uniquely fixed by the representation of the dynamical symmetry group. 
Every completely symmetric representation of SO(d+ 1) appears once and only once and 
corresponds to the multiplet with energy Eg. The angular momentum content of this mul- 
tiplet is determined by the branching rule fll9|). In three dimensions this expresses just the 
well-known fact that for each value of E n= £ + i the orbital angular momentum can vary from 
to n - 1 = I. 



III. N = 2 SUPERSYMMETRIC QUANTUM MECHANICS 

We wish to further generalize our results to the supersymmetric hydrogen atom in d dimen- 
sions. For that purpose we need a supersymmetric extension of <i-dimensional Schrodinger 
operators and in particular of the operator in (f|). Such supersymmetric Hamiltonians can 
be written as 



H = {Q, Q f } = #t with q2 = gta = , (20) 
where the supercharge Q and its adjoint 0} anticommute with a self-adjoint idempotent 
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operator T. The subspace on which T = 1 is called the bosonic sector and its complement 
the fermionic sector. Hence, Q transforms bosons into fermions and vice versa. From ( |20| ) 
one sees at once that the supercharge commutes with the supersymmetric Hamiltonian, 

[Q, #1 = 0, 

i.e. generates a supersymmetry of the system. The simplest models exhibiting the structure 
([20l) are 2 x 2-dimensional matrix Schrodinger operators in one dimension. Such models 



were first studied by Nicolai and Witten flHj, 11, |T2 



Supersymmetric Hamiltonians in higher dimensions have been introduced previously by 



several authors [y| |Tj]. Here we briefly present the construction used in this paper. We 



introduce a set of fermionic creation and annihilation operators, 

R,V>J} = <U, {ipaM = tyll>t} = Q > a, 6 = 1,..., <2 (21) 
and the Fock space with vacuum |0) which is annihilated by all operators ip a . This space 
splits into sub-spaces, 

C = C © d © . . . © C d , dimC p = f~J , dimC = 2 d , 

labeled by their 'fermion number' 

d 

N\ Cv = V 1 , where = ^ iplip a . 

a=l 

As basis in C p we may choose 

\a x ...a p ) ^ipl^.-ipljO) , ai < a 2 < . . . < a p . (22) 
Along with C the Hilbert space of all square integrable wave functions decomposes as 

H = H ®H 1 ® . . .®H d , where N\ Hp =pl. (23) 
An arbitrary wave function in 7i p has the expansion 

* = /ai...a p (z)K • • • a p) > f ai ...a p totally antisymmetric. (24) 
An explicit realization of the creation and annihilation operators can be given in terms of 
the hermitian 7-matrices in 2d Euclidean dimensions: ip a = \(j a ~ ijd+a)- 
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The supercharge and its adjoint 5 are defined via 



Q = e x Q e x = i^ ^ a (d a + d a x) , with Q = iip a d a , 

gt = e xgt e -x = J ^ ( a _ 9aX ) ? with gt = ^ . (25) 

a 

At this point the real superpotential x{ x u ■ ■ ■ i x d) remains unspecified. From (^l|) it follows 
at once that the free supercharge Qo is nilpotent and since Q and Qo are related by a 
similarity transformation the same holds true for Q. The supercharge Q only contains 
annihilation operators and hence decreases the fermion number by one. Its adjoint 
increases it by one, 

[N, Q] = -Q and [N, Q f ] = Q f . 

The supersymmetric Hamiltonian defined in (|20f) is a 2 d x 2 d -dimensional matrix Schrodinger 
operator and takes the following form (cf. also |U| |14}| ) 

iy = {-A + (Vx,Vx) + A X }l 2 <i-2^ViXa6^, Xab = -faj^ b - (26) 

a,fe=l 

We use brackets to indicate contraction of indices as (A,B) = ^2 a A a B a . Contrary to the 
supercharge and its adjoint the supersymmetric Hamiltonian H commutes with the number 
operator N and hence leaves each subspace Ti p in the decomposition (p3|) invariant, 

H : Tip — > Tip . 

On the subspace Ti p the supersymmetric Hamiltonian is still a matrix Schrodinger operator, 

H\ Hp = -Al + V^, tTl =Q- 

Only in the extreme sectors Tio and Ti p do we get ordinary Schrodinger operators acting on 
one-component wave functions. With 

V£V&|0) = and i>li> b \l2...d)=5 ab \l2...d) 

the corresponding potentials take the form 

^(o) = (V X ,Vx) + A X and = (V X , V*) - A X . 



5 The hermitian linear combinations Qi = Q + and Q2 = i(Q — Q') satisfy the standard TV = 2 
supersymmetry algebra = 2HSij. 
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More generally, for an arbitrary state \I* = f ai ...ap\cLi • • ■ ct>p) G H p the Hamiltonian acts as 
follows: 

p 

(a x ... a p \HV) = (-A + V^)f ai ... ap + 2 ]T (-)' Xai6 /^..^...^ . 



6,i=l 



(27) 



The nilpotent supercharges give rise to the following Hodge-type decomposition of the 
Hilbert space, 

H = QH © Q ] H © Ker H , (28) 

where the finite dimensional subspace Ker if is spanned by the zero-modes of H. Indeed, 
on the orthogonal complement of Ker H we may invert H and write 



which proves (p8"|). The supercharge Q maps every energy-eigenstate in Q^Ti PI H p with 
positive energy into an eigenstate in QH fl H p -\ with the same energy. Its adjoint maps 
eigenstates in QHC\H P into those in Q}H fl H p +i with the same energy. With the exception 
of the zero-energy states there is an exact 'pairing between the eigenstates and energies in 
the bosonic and in the fermionic sector as depicted below. 



Ho - — = — Hi - — = — H 2 



E 



Q 



QH 



Q 



w 



Hd-2 ' ^ ' Hd-i ' „ Hd 



QH 



Q^n 



QH 



w 



Q 















































QH 


Q*H 


QH 


Q^H 



Q 



The supersymmetric systems with supercharges fl25|) admit a generalized Poincare duality 
relating Hd- P with H p . This is seen as follows: instead of the vacuum |0) G Ho, which is 
annihilated by all i/; a , we take ||0) = ijr x • ■ ■ if> d \0) G Hd, which is annihilated by all ip\. As 



11 



basis in Cd- P we choose 

||ai . . .Op) = Va a • • •V'opllO) > 
such that an arbitrary wave function in 7id- p has the expansion 

* = /oi...oj,||ai • • • Op), /oi... ap totally antisymmetric. (29) 

The supersymmetric Hamiltonian acts on such a state as follows, 

v 

(O! . ..Op|| ff *> = (" A + " -^1 //„:....,..,,, ■ 

6,i=l 

A comparison with (^7|) yields the important duality relation 



IP 



H 



-x 



V 



(30) 



which states, that the Hamiltonians in 7i p and Hd- P are equivalent, up to a sign-change 
of the superpotential. Similarly on finds, that the action of Q on the state ( 22 ) in H p is 
identical to the action of on the state (|2"9"|) in TCd- p , up to a sign-change of \- 



IV. THE SUPERSYMMETRIC HYDROGEN ATOM AND ITS SYMMETRIES 

First we consider d-dimensional supersymmetric systems with spherically symmetric super- 
potentials and derive the conserved angular momentum J a b. The total angular momentum 
J a b is the sum of two terms: the orbital part L a b and the internal part S a b which transforms 
the components of a wave function. For a superpotential x{ r ) t ne supercharges simplify to 



Q = iipa(d a + x a f) and Q j = iipl(d a - x a f) , where f = — . (31) 



They should be scalars with respect to the rotation group SO(d). However, since d a is 
a vector and tp a a scalar with respect to L a b, we need to supplement orbital rotations by 
internal ones such that ijj a becomes a vector as well. One easily verifies that 

[Sdb, if>c] = K$acA ~ <W0a) , where S ab = -r(iplipb ~ rfV'a) , (32) 

and that the hermitian S a b satisfy the same commutation relations (|5]) as the orbital angular 
momentum. It follows at once from fl32| ) that the supercharge ( [31]) commutes with the total 
angular momenta 
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Jab — Lab + S a b , (33) 

since Q only contains scalar products of vectors operators. 

Next we prove that there exists a supersymmetric generalization of the Laplace-Runge-Lenz 
vector in d dimensions if the potential is ~ 1/r. For such a potential the supercharges (|31"|) 
do commute with the supersymmetric generalization of (^), 

C a = JabPb + PbJab + X a f(r)A (34) 

with a suitable operator A. One can show that / must be the function in (|3T| ) for C a to 
commute with Q. This function, along with A are fixed as follows: 

First, for C a to be a vector, the operator A must be a scalar under rotations induced by 
the J a b. Second, in the zero particle sector Ho the vector C a must coincide with C a in (||. 
Third, A should commute with the particle number operator since the J a b do commute. 
A general ansatz for A subject to these three conditions reads 

A = al- pN -j&S , S = x a ij a , x a = y } (35) 

with constants a, (3 and 7 which ought to be determined. Clearly, C a is a vector operator, 
such that 

[Jab, Cc] = K^acCfc — 8bcC a ) 

holds true. Let us now calculate the commutators between the C a and the supercharge. 



They should vanish for a suitable chosen function / in (34j) . We obtain 



(36) 



[C a ,Q] = 2{f^ b + f'Sx b }Jab + (3fx a Qo + ifx a {((3 + 1 )rf + 1 d r }S 

+ i{M a + fSXa} (1 - d ~ A) + i 1Xa f (d~N-l)S . 

The terms containing derivatives are 

fx a ((3 - 2)Q + 2ir(/ + rf)Sd a + i( 7 / - 2rf)Sx a d r . 

They cancel if 

/ = or x = -Ar and (3 = -7 = 2 (37) 

hold true. With this choice all but the terms proportional to S ao in the first line on the right 
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in (|36|) cancel, and we remain with 

iC a ,Q] = ^(a-d+l)(4> a -x a S) 

which vanishes for a = d—1. Hence, the supersymmetric extension of the conserved Laplace- 
Runge-Lenz vector reads 

C a = JabPb + PbJab - \x a A with A = (d - 1 ) 1 - 2N + , (38) 
and this is the main result of this section. The choice x = f° r the superpotential in 
(p5|) and (^) leads to the following supersymmetric extension of the Coulomb Hamiltonian 

tf = -A + A 2 -M. (39) 

Restricted to the zero particle sector this is the Hamiltonian of the hydrogen atom 6 and 
restricted to the d particle sector it corresponds to the electron-antiproton scattering system. 
The corresponding supercharge and its adjoint take the simple form 

Q = Qo- iXS and Q f = Q\ + i\S j , 

where the free supercharge has been defined in ( p5|) and the operator S = x a tp a has already 
been used in ( |35|) . 



V. ALGEBRAIC DETERMINATION OF THE SPECTRUM 

We proceed as we did in the purely bosonic case and calculate the commutator of two C a : 

[C a , C b ] = -AU ab (-A - ^) © -AU ab (H - A 2 ) . 

Up to the shift in H and the replacement L ab J ab this is the same relation as (^). The 
total angular momenta J ab and the vector operator 

Ka = v/ 4(A^- H) ( 4 °) 
form a SO(d + 1) algebra on the subspace of bound states (E < A 2 ). Finally we should 



6 We have to identify 77 = X(d — 1). The additional shift A 2 makes the lowest eigenvalue of this operator to 
be equal to zero. 
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calculate C a C a . If we can express this scalar operator in terms of H and operators that 
commute with H, similarly as we did in (|T0D, then we may solve for H. However, one soon 
realizes that this is impossible by only using the operators 1, N, J a bJ a b and H. However, we 
can express C a C a in terms of 1,N, J a bJ a b and the two operators QQ^ and Q}Q as follows: 

C a C a = 4(A 2 - H)K a K a = -2\ 2 J ab J ah + (2 J ab J ab + (d-2N- l) 2 ) QQ* 

(41) 

+ (2 J ab J ab + {d-2N + If) Q^Q , 

and this relation is sufficient to obtain the spectrum of the supersymmetric hydrogen atom. 
Each of the three subspaces in the decomposition ( p8|) is left invariant by the Hamiltonian 
and we may consider H on each subspace separately. Since Q 2 = we find H\qu = QQ^ 
and H\q^ h = Q^Q and can solve ([41]) for H in each of these subspaces 

H\ QH -QQ -X - {d _ m _ 1)2 - 4C(2) , (42) 
HI OtO-A 2 (rf-2iV+l) 2 A 2 

H\ Qfn -QQ-X - {d _ m + 1)2 - 4C{2) , (43) 

where C( 2 ) is the second-order Casimir of the dynamical symmetry group SO(d+l), 

C(2) = ^JabJab = -^JabJab + K a K a . 

All zero modes of H are annihilated by both Q and Q\ and according to ([11]) the second- 
order Casimir must vanish on these modes, such that 

C m\KeiH = • 

We conclude that every normalizable zero mode $ of H must transform trivially under the 
dynamical symmetry group, Jab^ = 0. 

To obtain the bound state energies we need to determine those irreducible representations 
of the dynamical symmetry group which are realized in 7i and the corresponding values of 
the second-order Casimir operator. The degeneracy of an energy level is then equal to the 
dimension of the corresponding representation. 

We use the abbreviation T>^ to denote multiplets of the orthogonal groups corresponding to 
Young tableaux of the form 

T 
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since in the following only those representations will appear. Let us assume, that each 
component function f ai ... ap (x), entering the state \& G H p in ([?§), is a harmonic polynomial 
of degree £, that is 

fai...a p (%) ^ ^ fa\...a p b\...b l %b\%b2 ' ' ' ^b t 

61,62,— 

with / symmetric in the b— indices and of zero trace in each pair of them. Since f ai ...a p 
is completely antisymmetric in the a— indices, these objects transform according to the 
completely antisymmetric representation 

r 



of the rotation group SO(d) generated by the S a b- On the other hand, each homogeneous 
polynomial f ai ,.. a (x) transforms according to the completely symmetric representations 



1 






i 



of the rotation group generated by the L &. It follows, that the wave function \1> G H p 
transforms according to the tensor-product representations 

Vl ® V[ = V l p _ x © V^ 1 © V^ 1 © V' +l . (44) 

Recall, that p is the fermion number and i the order of the homogeneous polynomials. 
For ip G 7i or ip G TCa the first factor on the left hand side in (pE4j) becomes the trivial 
representation and we only obtain the fully symmetric representations T>f on the right hand 
side, in agreement with our earlier results in the purely bosonic case. In the sectors TC\ and 
H-d-i the first representation on the right hand side of ( pB| ) is absent. For linear functions 
with t = 1 the second representation on the right is missing. Finally, when using the results 
(fH|), one should keep in mind that the representations with p and d — p of SO(d) are 
equivalent, ~ , and that for even dimensions the representations X>1 2 are reducible 
and contain one self dual and one anti-self dual multiplet. All representations of the rotation 
group SO(d) appearing in the tensor product fliif) should group together into multiplets of 
the dynamical symmetry group SO(d+l). To continue we need the following rules for the 
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branching of SO(d+l)- into SO ((^-representations 7 , 



Vi 



vi 



V 



p-1 



v e-i 



50(d) 



• (45) 



" SO(d+l) 

Now it is not difficult to see, that in the sector TC P all polynomials up to order £ appear in 
the branching of only two SOid + l)-multiplets up to two extreme representations, 



vt@v 



{vl ® (i e vl © . . . © x>fl - v^ + x>< } 



50(d) 



50(d+l) 

Of course, for p = 1 the last representation of the rotation group is absent. There is one 
notable exception to these branching rules for even d: in the middle sector 7i n =d/2 the correct 
branching rule reads 



Vl © VI 



{Vl ® (1 © Vl © . • • © vi) - v e n +l - V l n - V l n _ x } 



SO(d) 



SO(d+l) 

We summarize our results: In odd dimensions d = 2n + l the following representations of 
SO(2n+2) containing bound states arise in the various sectors of T~C for i > 1: 

Q f Q 



7~lc 



V[ 



v[ 
v" 



vi 



V 

Vi 



n-1 



V 1 



In all sectors but Ho we have £ G N. In the subspaces Ho we have £ G No and £ = 
corresponds to the trivial representation. The sectors H p>n support no bound states (see 
below) and therefore it suffices to consider the sectors with N < n. 

In even dimensions d = 2n the following representations of the dynamical symmetry group 
S0{2n+l) arise for I > 1: 



7 They can be obtained from jTa or by using the program LiE. 
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Hq 



Hi 



Q ] Q 



H n -i 



H r 



V' 



u n~l 
V 1 5= 



The values of the quadratic Casimir operator of SO(d + 1) entering the formulae ( |42]r43| ) for 
the energies are 

C (2) = d(£+p-l) + £(£ - 1) - p(p - 1) . (46) 

The dimensions of the representations T> l p are given in the appendix. In odd dimensions one 
uses the formula flAlj) with 2n replaced by 2n = d+1 and in even dimensions the formula 
( |A17|) with d+1 = 2n+l. In additions one must set 



(£ u £ 2 ,...,£ n ) = (£,!,..., 1,0..., 0) 

in these formulae. 

Now we are ready to determine all eigenvalues of the supersymmetric Hamiltonian ( p9|) with 
the help of the results ([4^,[|!|) and fl46|) as follows: 

In Ho only the symmetric representations T>( of S f O(c?+ 1) are realized. Since in addition 
Q\n — we obtain the following eigenvalues for H in ([42]) 

d- 1 



l+2^ y 

The index at the energy E indicates the O-particle sector. Since the supercharges commute 
with the dynamical symmetry group the multiplet T>[ is paired with the same multiplet in 
H\. The eigenfunctions in Hi are obtained by acting with on those in Ho- According to 
our previous results (see the figure above) there exists the additional multiplet T>^ in H\. 
This is obtained by acting with Q on the same representation in the two-particle sector. 
Hence H = QQ^ on this second multiplet and we obtain 

d-l 



Ei (V() 
Ei (Di) 



A 2 - 
A 2 - 



d - 1 + 2£ 

d-3 
d-l + 2£ 



e n 



e n 
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Note that t = does not occur in 7i\. In Ho the states with i = have vanishing energy 
and hence are annihilated by . 

Now one continues with 7i 2 and further on to H3 etc. One only needs the formulae 



rip 



d + 1 - 2p 
d-l + 2£ 
d-l-2p 
d-l + 2£ 



A 



g n 



G N 



We shall determine the corresponding eigenfunctions in the following section. 



VI. EIGENSTATES OF THE SUPERSYMMETRIC HYDROGEN ATOM 

So far we have not considered which highest weight states of the dynamical symmetry group 
are normalizable. Now we explicitly construct these states in all subspaces H p C H. In 
the previous section we have seen that for any £ > 1 there are one or two irreducible 
representations of SO(d + 1), namely 

V l p c (H p n Q ] H) and C (H p n QH) . (47) 

It suffices to construct the highest weight states ty p (V£ +1 ) of the latter multiplets. The 
highest weight states ty p (U£) of the first set of multiplets in ( fj?j) are then just their super- 
partners, 

Actually we only need to calculate the highest weight states ^ P (T> p+1 ) for p < d/2 because 
of the duality relation 

(p,A) < — > (d-p,-X) , 
which follows from (|30|) and (j37[). 

Observe that for any normalizable if-eigenstate ^ G QH the transformed state is 
normalizable, as can be seen from 

(Qttf , Q f ^) = QQ^) = (*, fftf) = , . 

Without calculating the highest weight states we can argue in which sectors bound states 
cannot exist. For that purpose we consider the Hamiltonian (|39|). It is easy to see that the 
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hermitian operator S^S, where S has been defined in (|35|), is an orthogonal projector, and 
hence has eigenvalues and 1. It follows at once that for p > d/2 the operator A in (|38| ) is 
negative and hence H > A 2 . We conclude that H has no bound states in the sectors TL p> d/2- 
On the particular sector 7i n the operator A has both positive and negative eigenvalues. We 
expect that in this sector only one of the two representations (for each £) of the dynamical 
symmetry group contains bound states. After these general considerations we proceed with 
computing the highest weight states ty p (Vp +1 ) in the subspace TC P fl QH. Again we proceed 
differently in even- and odd-dimensional spaces. 

Even-dimensional spaces: We use the complex coordinates Zi,...,z n in R d=2n and the 
creation/annihilation operators <p\, . . . , <$ w <pi, . . . , <fi n introduced in appendix |A 1| , together 
with the complex conjugated objects. Since the dynamical symmetry group is SO(2n+l), we 
should take the Cartan operators if; and raising operators E,i from appendix ( |A 2|) . However, 
we must remember that the last row and last column of (Jab) contain the components of 
the generalized Laplace-Runge-Lenz vector. As a consequence the last step operator E n in 
( |A16|) is to be replaced by 



E n = 4|(^-i + iK d ) ~ -4(Cd-i + iC d ) 

Xz ( 48 ) 
= -2z n A + (2rd r + d- 1)4 - 20t(M + fadj + 2(0^ + 4Wn ~ , 

Since the simple roots of the rotational subgroup SO(d) are positive roots of the dynamical 
symmetry group, a highest weight state of SO(d+l) is automatically a highest weight state 
of SO(d), similar as in the purely bosonic case. Since the two groups share the same Cartan 
generators the highest weight state ^ p {Vp +1 ) of SO(d+l) must also be a highest weight 
state of the multiplet T>^ +1 of SO(d). From the branching rule ( ^5|) and the tensor products 
([44]) it follows, that this highest weight state must be the state y a (£,P+ 1) given in (|A13|) . 
Hence we are lead to the ansatz 

%(V* +1 ) = f(r)y a (£,p + l). 

It remains to determine the radial function f(r) such that is annihilated by E n . With 



E n from one finds the following equation for the radial function /: 

(d - 1 + 2£)f + X(d - 1 - 2p)f = 
such that the relevant highest weight states in the p-particle take the form 
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d — 1 — 2p 

^p( V p+i) = ex P (-7*0 3>a(4 P + 1) with 7 ^ = - _ - - ^ A . (49) 
As A is assumed positive, these are bound states for p < n. 

Odd-dimensional spaces: For odd dimensions d = 2n + 1 the rank of the dynamical 
symmetry group SO(2n + 2) exceeds the rank of the rotation group SO(d) by one. As 
in the purely bosonic case we combine the first 2n coordinates, creation and annihilation 
operators to complex ones (cf. appendix A~2 ). Since the rank of the dynamical symmetry 



group is even, we should take the Cartan generators from appendix |A 1| with n replaced by 
n + 1. Since (Jab) contains the Laplace-Runge-Lenz vector the explicit realization of these 
operators differs from the one in this appendix. More precisely, the first n Cartan generators 
are those in ( ]A4jj , but the last one H n+ i is K d ~ C d (cf. appendix |TJ), where 

C d = -2x d A + (2rd r + d- l)3 d - 2^J(M + 4>A) + 2(^3 + 4>\ B i)^d - \x d A . (50) 
The raising operators are the n — 1 operators in ( |A^) plus the two operators 8 



E n = -r(z n d Xd - x d d Zn + 4> ] n i) d - ^\(j) n ) and E n+1 = -J=(C d _ 2 + iC d _i) . ( 51 ) 

The last operator coincides with E n in (0). By using similar arguments as in even dimen- 
sions we are lead to the following ansatz 

%(v* +1 ) = f(r)y a (e, P + i) 

for the highest weight state of T>^ +1 C Tip. This function is annihilated by all Ei< n . The 
condition E n+ i^ p = yields the same differential equation for the radial function / as before 
and we obtain 

d — 1 — 2p 

*i>fti) = exp (-7V) y a {£, P+l) with llv = d _ 1 + 2£ A • (52) 

For positive A these states are normalizable for all p < n. It is easy to see that the last 
Cartan generator ~ C d annihilates this state and this shows that it has the correct weight. 
The remaining highest weight states: We have argued that the highest weight state 
^/ p+ i(Vp +1 ) C Hp+i is the superpartner of ^ p {V^ +l ) in (p2^). A simple calculation yields 



which arc independent combinations of the last two raising operators, see appendix H. 
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(£+p)y s (£,p+l))exv(-y lp r) 



(53) 



for this state and shows that it is a linear combination of the two highest weight states y s 
and y a of SO(d) given in formulae ( |A12| , [A13| ). These states lead to additional bound-state 
multiplets in the sectors TC P with p = 1, . . . , n. 



VII. THE SUPERSYMMETRIC HYDROGEN ATOM IN DIMENSIONS < 4 



In this section we apply the general results of the previous three sections to supersymmetric 
systems in low dimensions. The two-dimensional case mainly serves as illustration of the 
method. It may be worth noting that such systems admit a supersymmetric Laplace-Runge- 
Lenz vector, contrary to what has been claimed in the literature |TB[ . The three-dimensional 
quantum system is of course the most interesting extension of the ordinary hydrogen atom. 
We have included the four-dimensional supersymmetric system since it already shows very 
nicely which additional structures arise when one goes beyond three dimensions. 
Two dimensions: The Hilbert space of the supersymmetric hydrogen atom in two space 
dimensions contains three sectors, 



n = n © Hi © n 2 , 

and an arbitrary wave function has the expansion (in the basis (|22|)) 
* = /o|0) + (/i|l) + / 2 |2» + /i 2 |12> ~ / = (/ , fi, h, fi2) T - 



The Hamiltonian (|39D acts on the component functions in / as follows 

/-l 



/ — >Hf with H = -A + A 2 + 



A 



S ab - 2x a x b 

V o o 1/ 

Clearly, for A > there are no bound states in the two-particle subspace, in accordance with 
our general result below eq. (|47|) . Only the multiplets 



V^°CH and ^ >u C^i 
of the dynamical symmetry group SO (3) contain normalizable states. 
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We introduce the complex coordinate z and the complex creation/annihilation operator (see 
appendix |A 1| ) in terms of which the highest weight state (^) read 

A 



M^D = ex P (- leo r)y a (£, 1) = exp(- 7 «,r) z"\0) 
and its superpartner ( |53"D 



l + 2£ 



*i(X>/) = Q^o(V() = i{(\- leo )Sly a (£, 1) + ey.(£, 1)J e X p(- 7 , r) . (54) 
The energy of the £+1 states in each of the two corresponding SO (3)-multiplets is 

Et = A 2 - 7 , 2 o • 

There is exactly one zero-energy ground state in the zero-particle sector and this state has 
I = 0. It is annihilated by the adjoint supercharge Q\ as can be seen from ([54]) . The 
spectrum of the supersymmetric system is depicted in the following figure. 



7^o 



Hi 



E/X 2 
1 



Eigenvalues of H in d = 2 dimensions. 
Three dimensions: For this most relevant system, the Hilbert space contains four sectors, 

H = H © Hi © H 2 © n 3 , 
and a wave function has the expansion 

* = /o|0> + (/i|l> + M2> + /s|3» + ... • 

Since here we are only interested in bound states it suffices to consider H on the subsectors 
H.q and Tii. The Hamiltonian in T~C belongs to the ordinary hydrogen atom, 

2A 

r 



H (0) 



-A + A 2 
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and it acts on a state \I/ G Tii with component functions (a\H) = f a as follows, 



(a\HV) 



2A 

A + \ 2 )f a X a X b fl 



b ■ 



We take the coordinates (z, X3) and the creation operators (0^, ipl) (see appendix |A2j ). 
The highest weight state in the multiplet T>[ C H.q of the dynamical symmetry group SO (4) 
has the form 

V (V() = exp(- lm r)y a (£, 1) E Ho , 7£0 
and this state is mapped into the partner state 



A 

IT 



*i(X>fl = Q^oOP/) = if (A - 7^3^, 1) + ey,{£, 1)) exp 



All states in the two corresponding 5'0(4)-multiplets share the same energy 



Ei — X — Ym 

and both multiplets contain (£ + 1) 2 states. This is the well-known spectrum of the hydrogen 
atom. The normalizable zero-mode has I = and resides in the zero-particle sector. It is 
just the ground state of the hydrogen atom. 

There are no further bound states, since the other states are paired with wave functions in 
Ti.2 and thus cannot be normalizable. The spectrum of the supersymmetric hydrogen atom 
in three dimensions is shown below. 



7~lc 



n. 



E/X 2 

1 



Eigenvalues of H in d = 3 dimensions. 



Four dimensions: The Hilbert space splits into five subsectors, 
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In the zero-particle sector the Hamiltonian takes the form, 

3A 



# (o) = -A + A 



2 



r 



and in the one-particle sector it acts on a state \1/ with (a\^) = f a as follows: 

(a\m) = (-A + A 2 - -) f a - —x a x b f b ■ 
\ r J r 

In Ti.2 we are left with a 6 x 6-matrix Schrodinger operator: For a two-particle state ^ with 
component function {ab\^/) = f a b we obtain 

{ab\HV) = ( -A + A 2 + - ) f ab - —(x b f ac - x a f bc )x c . 



We introduce complex coordinates Zi,z 2 and creation operators as m appendix |A 1 

For all £ G N we find the highest weight states 

, s i , „i . . , 3A 



# (£>f) = exp(- 7a) r) z[ |0) with j eo 



2£ + 3 

in the zero-particle sector, together with their superpartners in Tix, 



^ 1 (V() = i((A - 7io)&y a (£, 1) + £y s (£, 1) j exp(- 7£0 r) , 
which exists for £ > 0. All states in these two multiplets have the same energy 

^ = A 2 - 7 , 2 , 

and the number of states in each multiplet is (£ + + 2)(2£ + 3)/6. Now there is an 
additional representation £>| C Hi with highest weight state 

^(D-f) = exp(- lei r)y a (£,2) 

= exp(-7ar) {4>\z x - 4>\z 2 ) z[~ l \$) , 7a 



3 + 2^ 

The energy of this state is 

^ = A 2 -7a, 

and the multiplet contains £(£ + 3)(2£ + 3)/4 members. Again there is a supersymmetric 
partner multiplet with the same energy and degeneracy. The remaining highest weight state 
of T>1 ~ T>2 C 7^2 is paired to a state in and therefore not normalizable. Thus we find 
the spectrum as depicted in the following figure. 
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Eigenvalues of H in d = 4 dimensions. 

VIII. SUMMARY AND DISCUSSION 

We have succeeded in extending the celebrated results of Pauli, Fock, Bargmann and others 
in two directions: to higher dimensions and to the M = 2 supersymmetric extension of the 
hydrogen atom. First we constructed a generalized angular momentum and an extended 
Laplace-Runge-Lenz vector which could be combined to generators of the dynamical sym- 
metry group SO(d+l) in d dimensions. Then we related the quadratic Casimir operator 
of this group to the particle number N, QQ^ and QQ^ . This way we calculated the bound 
state spectrum of the supersymmetric hydrogen atom in arbitrary dimensions by algebraic 
means. We have determined all bound-state multiplets of the dynamical symmetry group 
and calculated the highest weight state in each of them. 

Itzykson and Bander || distinguished between the infinitesimal and the global method to 
solve the Coulomb problem. The former is based on the Laplace-Runge-Lenz vector and is 
the method used in this paper. In the second method one performs a stereographic projection 
of the (i-dimensional momentum space to the unit sphere in d+1 dimensions which in turn 
implies a SO(d + l) symmetry group. It would be interesting to perform a similar global 
construction for the supersymmetric systems introduced in this paper. 
We have not explained why every multiplet of the dynamical symmetry group appears four 
times. Furthermore there is a new 'accidental' degeneracy: in higher dimensions some 
eigenvalues of the Hamiltonian appear in many different particle-number sectors. It may 
very well be, that the algebraic structures discussed in the present work have a more natural 
setting in the language of super algebras. We have not investigated this question. Anonther 
interesting question is whether the dynamical symmetries considered in this paper are related 
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to fermionic Killing- Yano supercharges |21| . This problem needs further investigations. 
Finally, we see no obstacle against extending our method to the scattering problem of the 
super symmetric hydrogen atom, for which the dynamical symmetry group SO(l,d) is non- 
compact. 
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Appendix A: REPRESENTATIONS OF ROTATION GROUPS 

In this appendix we collect the group theoretical facts needed in the main body of the paper 
(cf. |15|, [17], [18], for a more detailed discussion of these issues). We shall construct the 
relevant irreducible representations of the total angular momentum operators 

Jab = t [ Xa g^' b ~ Xb Q^~) + Y (fa^b ~ , a,b=l,...,d 

satisfying the so(d) commutation relations 

[Jab, Jed] = i($acJbd + hdJac — 5 a dJbc ~ $b C Jad) 

on wave functions in 



Hp = L 2 (R d ) x with p = 0,...,d . 

The fermionic operators tp a have been introduced earlier in section [TTT| It is convenient to 
use the Cartan-Weyl basis consisting of generators Hi in the Cartan subalgebra and one 
raising and one lowering operator E a and E_ a for every positive root a, 

[H h E a ] = aiE a and [E a , E- a ] = a ■ H with E^ a = E ] a . (Al) 

Because of their different properties we do this separately for the groups D n ~ SO(2n) and 
B n ~ 50(271+1). 
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1. Total angular momentum for the SO(d = 2n) groups 

To proceed it is very convenient to introduce complex coordinates in R 2n 



If . v If . , 

Zi — —7={X2i-l + lX2i) , Zi — —=\X2i-l — lX2i) , 

1 - 1 (A2) 

di = -j={d X3l _i - id X2i ) , di = —(d^^ + id xai ) , i = 1, . . . n 

and similarly two sets of complex creation- and annihilation operators 

1 - 1 (A3) 

4>i = -j={ip2i-i - ifai) , 4>i = ~^2^ 2i - 1 + x ^ 2 *) ' 1 = ' ' • n ■ 

The only non-vanishing anticommutators are 

{<t>i, 0j} = = &ij ■ 

The generators in the Cartan subalgebra take the simple form 

Hi = J 2 i-i,2i = Zidi - Zidi + - 4>\4>i > i = l,...,n (A4) 
and there are two types of raising operators: 

E a = ^(J2i-l,2i-l + J2i,2j - i^2i-l,2j + i^2i,2j-l) With TOOt a = Ci - Cj , 

E a = -(J2t-i,2i-i - hi,2j + iJ2i-i,2j + iJ2i,2j-i) with root a = ei + ej , 
where i < j is assumed. In terms of the complex coordinates/operators they read 

E a = — [zidj - zjdi + 4>\(j)j - 4>)<Pi) with root a = d - ej , 



±J a = — yZidj — Zjdi + <Pi<Pj — (PjVi) with root a = e, + ej . 

The corresponding lowering operators are just the adjoint of the raising operators. The 
operators (Hi, E a , E- a ) satisfy the commutation relations ( |AT| ) with corresponding positive 
roots in ( |A5|) . The n simple roots are 

o.i = ei — ej + i , 1 < i < n and a n = e n _i + e n 

and the corresponding raising operators have the form 

Ei = \ (zid i+1 - z i+1 di + 4>\(f)i + i - 4>\ + i4>i) , a = ei- e i+1 , 1 < i < n (A6) 
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^ ^n— 1 "I - • 



(A7) 



E n - - (yZ n _id n — Z n d n -\ + 

With the help of the Weyl vector 

5 = ^^a = (n- l)ei + (n- 2)e 2 + . . . + e n -\ , 

where the sum extends over all positive roots in (|A5|), we may calculate the dimension of 
an arbitrary faithful representation of SO(2n). Such a representation is determined by its 
Young tableau containing at most n rows. The length ^ of row % is bigger or equal to that 
of row i + Hence, a Young tableau is given by n non- negative ordered integers 

h > h > ■ ■ ■ > 4-i > L 

and has the form 



P 



1 


2 








ii 


1 


2 






h 
















1 




£„ 







, p < n . 



Rows with length are not shown when one draws a Young-tableau. The corresponding 
representation p^i>->*« has the dimension 



dim (Z>*~^) = \\ 



+ £ s + 2n-r-s £ r -£ s + s-r 



2n — r — s 



s — r 



l<r<s<n 

For the second-order Casimir invariant of these representations one obtains the formula 

e {2) (V^- In ) = J2 l r(tr + 2n- 2r) . 

r 

In particular, for the completely symmetric representations 

v m...o _ v t 

these formulae simplify to 

C (2) (V() = £(£ + d-2) and dim(P/) = 
For the completely antisymmetric representations 



(A8) 



1 


2 






£ 



'+d-3 
£-2 



(A9) 



pi,i,..i = pi „ 
p 
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they simplify to 

C {2) {V$)=p{d-p) and dim(P p 1 ) = 
Simultaneous eigenstates of all n generators ifj in the Cartan subalgebra have the form 

i=i 

where mi,rhi e No and G {0, 1}. The vacuum |0) is annihilated by all particle lowering 
operators ?/> a or equivalently by all fa and 0j. The ifj-eigenvalues of these states are m — 
fhi+pi-Pi. 

Next we must construct the highest weight states which are annihilated by all raising oper- 
ators. Every such state determines an irreducible representation. The eigenvalues of Hi on 
a highest weight state is equal to the length ^ of the Young tableau corresponding to the 
irreducible representation determined by this weight. The d + 2 space-independent highest 
weight states are 

\p) = \pp) with Pi > ■ ■ ■ > Pn > p' n > ■ ■ ■ > p[ 
and ^(Pi + p'i) = p , 

There is an additional highest weight state in the p = n particle sector, that arises since in 
this sector we have selfdual and anti-selfdual configurations. It is given by 

Pl = ...= p n _ x = p' n = 1 , p n =p' 1 = ...= p ' n l = . 

Clearly, the particle number p uniquely specifies these state since the Pi,p'iS are ordered. 
These states define the completely antisymmetric representations 

Dp for p < n and ~ T^\ n - P for p > n . 

We used that a Young tableau, the first column of which has length n < p < 2n, gives rise to 
the same multiplet as the tableau with first column of length 2n—p < n. It the following one 
should replace T>^ by T>2 n _ p if p exceeds n. Also note that T>q ~ T>\ is the one-dimensional 
trivial representation. 

The highest weight states in the 0-particle sector are 

*f|0) (A10) 
and they give rise to the completely symmetric representations T>[ spanned by the harmonic 
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polynomials of order £. The relevant irreducible representation of SO(2n) in the p-particle 
sector is gotten by tensoring the antisymmetric representation D p x with a symmetric repre- 
sentation Df. We use 



Vl ® D[ = D" © D"- 1 © X? m © D, 



p+i 



(All) 



or in the language of Young tableaux, 



1-1 



1 2 



p — 1 



1 








1 


2 








V 


© 


V 





>p+l . 



Note that for p—1 and/or I = 1 there appear only three representations in this decomposi- 
tion. For p—1 the first representation and for i = 1 the second representation on the right 
hand side in ( [All ) are absent. Also note that for p = n the first and last representations 
are equivalent. The second to last representation Dp +1 on the right hand side has highest 
weight state 



y s (£+l,p)=z{\p) , (A12) 

as it is the product of the highest weight states of Dp and D[. To find the highest weight 
state of the other representations we observe that the operators 



rS = x a 4> a = Zi4>i + Zifa and rS 1 = x a ^\ = + Zi<j)\ , 



which have been introduced in (pop, commute with the total angular momentum and hence 
map highest weight states into highest weight states. Since 5* decreases and S> increases the 
particle number by one, we find the state 



y a {£,p + i) = rsy s {£,P+i) 



p+i 

i=l 



z i 4...ti...<t>l +1 z{- 1 \o) 



(A13) 



which is highest weight state of the last representation D^ +l in the decomposition ( [All ). 
The missing two highest weight states correspond to those representations in the tensor 
product of a symmetric and an antisymmetric representation which one obtains by taking 
the trace over two suitable indices. This operation is equivalent to acting with S'. Thus 

T.(e,p-i) = &y.(£,p-i) 
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is the highest weight state of T>^_ x in the decomposition ( |A11| ). For the remaining highest 
weight state we make the ansatz 

r a (i-i, P ) = (ss^ + atfs)y a (£-i,p) . 

As {S, S^} = 1 this state may have a component in the direction of y s (£ — l,p). However, 
for the choice a = —1 the highest weight state 

T a {i-i,p) = [s,s*]y a (i-i, P ), 

is orthogonal to y s (£—l,p). 

2. Total angular momentum for the SO(d = 2n + 1) groups 

The rotation group SO(2n+l) has the same rank as its subgroup SO(2n) and hence we may 
still use the Cartan generators (|A4|) , that is 

Hi = J 2 i-i,2i = - Zidi + 4>\(j)i - 4>\(j)i , i = 1, . . . , n . (A14) 

We use the complex coordinates and the complex creation- and annihilation operators 
(|X3|), supplemented by the last coordinate Xd and the last creation and annihilation operator 
ijj d and ipa- Clearly, the raising operators ( |A5| ) are still raising operators of so(2n + 1) with 
the same positive roots. But since 

dim(SO(2n + l)) = dim(SO(2n)) + 2n and rank (SO(2n + 1)) = rank (SO(2n)) 

there are n positive roots missing. These are 

E a = (J 2 i-i,d + U 2i ,d) = t (zid Xd - x d di + (f)\ip d - ijj d 4>i) , a = ei , 

where 1 < i < n. The first n — 1 simple roots are the same as in ( |A6|) , but the the last one 
is replaced by e n . Hence the raising operators corresponding to the simple roots read 

Ei = t (zid i+ i - z i+1 di + 4>\<p i+ i - 4>\ +1 4>i^ , a = ei- e i+1 , 1 < i < n , (A15) 

E n = t (z n d XA - x d d n + 4>l4> d - i) j d 4>n^ , a = e n . (A16) 

The Young tableaux are identical to those of SO(2n) and hence are characterized by n or- 
dered non-negative integers Ei, . . . ,£ n . The dimensions of the corresponding representations 
read 
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j . ,~. t . e\ ~fj2£ t + d-2t -i-r £ r + £ s + d-r-s£ r -£ s + s-r 

dim(v )=n d _ 2t n — — — (M7) 

t=l l<r<s<n 

and the formula for the second-order Casimir is the same as for the so(2n) algebra, 

C {2 )(V £l = J2 £ r(tr + d-2r) . 

r 

Also the rules for tensor products are identical to those of SO{2n). 

Since the simple roots are different, the highest weight states have a slightly different form. 
The simultaneous eigenstates of the n generators in the Cartan subalgebra read 

f( Xd ) J] z^zf 1 \pqf) , \pqp') = 4 P1 ■ ■ ■ € Pn ^\ q A A ■ ■ ■ ft^lo) , 

i 

where m^m, G Z and Pi,q,pi € {0, 1}. The d+1 constant highest weight states are 

\p) = \pqp') with p x > . . . >p n > q >p' n > . . . >p[ , 

where p = Yl(Pi~^P'i) + Q denotes the particle number. The highest weight of T> p +1 in the 
decomposition 

V x v ® V[ = V l v _ x © V^ 1 © V l p +l © V l p+l (A18) 
is again determined by the highest weight state 

y s (£ + i, P ) = 4 \ P ). 

As in even dimensions one may use the scalar operators 

rS = x a i) a = Zi4>i + Zi4>i + x d i) d and rS ] = x a ^\ = z^\ + z^\ + x d ip a 
to obtain the highest weight states 

y a (£, P + i) = rSy s (£,p+l) — v^ +1 , 
T s (£,p-1) = rfy a (£,p-l) — > V', 



r a (£-i, P ) = [s,s*]y a (t-i, P ) — v 

of the remaining irreducible representations in (A18). 



-l 
p > 
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Appendix B: ROTATION GROUPS VS. DYNAMICAL SYMMETRY GROUPS 



In the main body of the paper we have seen, that the total angular momentum J a b in fl33| ) 
together with K a in (HO) combine to generators of the dynamical symmetry group SO(d+l) 



J 



Jab 


K a 


-K b 






AB 



The rotational group with generators J a b discussed in the previous part of the appendix, 
must be embedded into the dynamical group, 

d = 2n: SO{2n) C SO(2n + 1) 

d = 2n+l : SO(2n + 1) c SO(2n + 2) . 

Even dimensions: The dynamical symmetry group has the same rank as the rotation 
group SO(2n) and we can repeat our construction in appendix |A2j , where we extended 
SO(2n) to SO(2n + l). Of course we should take into account that the components in the 
last column and last row of (Jab) are the components of K a . The Cartan generators are 
those in ( A14 ) and the first n — 1 raising operators are given in ( A15 ). But the last raising 
operator ( A.16 ) is of course replaced by 

1 



En 



V2 



(K^ + iK d ) , 



which is proportional to ^(C^-i + iC^). The latter has been given in (|48|). 
Odd dimensions: The rank of the dynamical symmetry group SO(2n+2) exceeds the rank 
of the rotation group SO(2n+l) by one. The Cartan generators are given by the n operators 
Hi in (|A4p, supplemented by H n+ i = K d ~ C d , where the explicit realization of C d is given 
in (pD|). The raising operators are the n — 1 operators E\ in ( |A15| ) plus the two operators 

1 



E n 



E' 



(J 



d-2,d 



+ K d _i - ii^ d „2 + Ud-iA) 



(Jd-2,d — Kd-i + iKd-2 + iJd-i,d) 



a 



a 



G71 Cn+l ; 



Highest weight states are annihilated by these two raising operators and it is convenient to 
use two (independent) combinations of these operators, namely the operators 



■±=(E a + E' a )=E n and (E a — E' a ) ~ -^=(C d _ 2 + iC d -i) 



E 



n+l 
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Their explicit forms can be found in (f48|) and (|51|). 
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